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In this work there is given an approximate solution of the problem 
of a two-dimensional steady supersonic stream of ideal gas, neglecting 
heat conduction, flowing around a thin wing with sharp edges at small 
angles of attack. (Determination of the law of distribution of pressure 
along the wing, lifting force and head resistance of the wing.) 


PART I 


The problem of the investigation of the mechanical action of a 
moving gas on an immovable wing appears as a special case of the some- 
what more general problem of the investigation of the mechanical action 
of a moving gas on an immovable fixed wall constraining the motion of 
the gas. In our own explanation we begin with the formulation of this 
last problem in which we confine ourselves only to the consideration of 
the steady two-dimensional forces of ideal gases not subject to the action 
of gravitational forces. In the plane of motion of the gas we shall 
arrange an immovable rectangular coordinate system in such a manner that 
it is situated as in figure 1. We introduce three functions v, p, anri p 
of the independent variables x and y defined, respectively, as the veloc- 
ity, density, and pressure. The vector functions v will be determined by 
a pair of scalar functions of the independent variables. For these func- 
tions we shall agree to take either the functions v x , Vy defined as the 

projections of the velocity of the axis x and y, respectively, or the 
functions v and P, defined, respectively, as the absolute value of the 
velocity and its angle with respect to the positive direction of the x-axis, 
measured in the counterclockwise sense. In what follows we limit ourselves 
to the consideration only of flows for which the function p satisfies the 
condition 


- — < B < — 

2^2 


( 1 ) 
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As is well known, the study of the gas motion under cons ider at ion 
leads to the investigation of the following system of differential 
equations 


^x 5v x t 1 dp 

dx ^ dy p dx 




0 


' Tx !^ + Vv ii + iap,o 

dx dy P dy 

> 

S(PV X ) 5(pVy) _ 
dx dy 


dxln^J 





( 2 ) 


Here k is the adiabatic exponent (for air k = 1. 405) • If the motion 
of the gas is constrained by an immovable frictionlesB fixed wall in the 
plane XOY, the gas will be adjacent to it along some curve. We shall 
call this curve the "contour K." 


Consider the unit vector t tangent to the. contour K directed in 
such a manner that its projection on the x-axis is positive. Denote 
by the angle which it makes with the x-axis. Clearly may be 

regarded as a function of tlje abscissa x of that point of the contour K 
associated with the vector t. We denote this function by Pjj(x) and 

assume that it is continuous. If the function p^(x) is prescribed and, 
moreover, the coordinates of any point of the contour K are given, the 
form and position of the contuur Is completely determined. We agree to 
take as origin the left edge of the contour K. Then the equation of 
this contour will have the form 


-Jr’ 




y = P tan p k (x)dx 
J 0 


( 3 ) 
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We can write this equation more briefly if we designate its right-hand 
side by y k (x) 

y = y k (x) (4) 

Since in the flows under consideration the direction of the velocity on 
the contour K must coincide with the vector t, the condition on the flow 
along an immovable fixed frictionless wall may be written in the following 
fashion 


P = Pk( x ) 


(5) 


at y = y k (x) . The condition ( 5 ) must be added to the system of equa- 
tions (2) as a qualifying boundary condition. Much work has been dedicated 
to the investigation of solutions of the system (2) subject to the con- 
dition ( 5 ) • Of these we are interested here only in those in which the 
flow is supersonic, i.e. , flows at every point of which the following 
condition 


v > a 


(6) 


is satisfied, where a is the local speed of sound 



The investigations contained in these works divide in two fundamental 
directions. The first direction is represented in works in which solu- 
tions of the problem are achieved with the help of numerical or graphical 
processes permitting the step-by-step calculation of a system of parti- 
cular values of the desired functions. (Works of Busemann, Kibelia, 
and Frankl.) The fundamental achievements of the methods represented by 
these works consist of the fact that by their use many actual practical 
problems may be solved quantitatively of which the solution by other 
methods would present great difficulties. In particular these methods 
solve thoroughly corner -nonpotential problems. The chief defect of 
these methods is that the solutions obtained are numerical so that it 
is Impossible to obtain a general qualitative estimate of the phenomena 
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under investigation. The second direction is represented by the works 
of Meyer, Ackeret, Prandtl, and Busemann, which are confined to a culti- 
vation of an exact theory of irrotational flows. The results are based 
on the fact that in the case where vorticity is absent the character- 
istic system of differential equations (2) admit of integrable combina- 
tions. This theory leads to series of approximate results of any desired 
accuracy, giving a complete qualitative and quantitative picture of the 
flow. Since our investigation is mostly connected with the theory of 
irrotational flows we give below a brief introduction to the fundamental 
methods and results of this theory. 

We introduce the stream function i|r defined by the following relations 


8x 




8y 


= -pv x I 


J 


(8) 


As is well known from equations (2), ( 7 ) , and (8) the following relations 
follow without difficulty 



(9) 


v£ 

2 


+ 



t 0 


( 10 ) 


Where 9, tQ denote quantities which display the flow once and for all 

as a function only of i)r. With the help of equations (2), (9) , and (10) 
it is easy to obtain the two equations 



8x dy 


= n 


( 11 ) 
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. dv. 


dv,. 


_ fa V . 


(a 2 - v x 2 )^ + (a 2 _ Vy2) pL . v x vJ — ^ + PL) - 0 


dx 


6y 


M — *• + — - 
\8x dy 


( 12 ) 


where 2 denotes a quantity defined as 


ft = p 


dt r 


d In 6 


dijr k(k - l) dijf 


( 13 ) 


Equations (ll) and (12) represent linear relationships between the first 
partial derivatives of the functions v x , Vy with respect to x and y. 

Since every flow under consideration is supersonic, the entire region 
of the flow may be covered by a pair of families of characteristics. The 
differential equations of these characteristics are obtained easily by 
the use of equations (ll) and (12) . For one family of characteristics, 
which we shall agree to call the first family, we obtain the equations 


dy = m-jdx 


m 


d(v cos p) + ra. 2 d(v sin 0) = ft 


(a 2 - v2cos 2 p)m]_ + v^sin p cos p 

v 2 cos2p _ a 2 


and for the other, which we shall agree to call the second family, we 
have the equations 


dy = mgdx 


( 16 ) 


/ . / V (a2 - v2cos2p)mp + v2sin p cos p , 

d(v cos p) + mqd(v sin p) = ft — dx ( 17 ) 

v2cos2p - a 2 


Here mq, mg denote the following expressions 


mq = 


-v 2 sin p cos p + a \^v 2 - a 2 


a^ - v^cos^p 


( 18 ) 
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m.2 = 


-v2sin p cos p - a \/v2 - a2 

a^ - v^cos^p 


( 19 ) 


We now consider a supersonic stream with constant hydrodynamical 
elements (i.e., functions y, p, p, p, a). We shall call this flow the 
undisturbed flow. The values of' the functions v, p, p, a in the undis- 
turbed stream will be denoted by w, Pq, Pq, ag respectively, and 

the ratio w/ao by M. Since the stream under - consideration is super- 
sonic, M > 1. We shall choose the direction of the velocity of the 
undisturbed stream to correspond to the direction of the x-axis. 

We assume that the undisturbed stream strikes an immovable, fixed, 
frictionless wall (contour K) , inclined in such a manner that in flowing 
around this wall the stream never detaches from it and -remains super- 
sonic everywhere. We mdy distinguish two cases of flows of this type. 

Case I .- The contour K is situated in such a manner that the 
condition 


Pk(0) £ 0 (2°) 


is fulfilled. In this case, as is well known, there appears a curve of 
weak discontinuity OC (figs. 2(a), and 2(b)) proceeding from the origin 
and dividing the entire flow in two parts. On one side of the curve of 
weak discontinuity OC extends the region containing the undisturbed 
stream and on the other the region of flow around the wall. In the region 
of flow around the fixed frictionless wall the hydrodynamical elements of 
the stream, generally speaking, are not constant but vary. In what fol- 
lows we shall call this part of the stream the disturbed stream. In the 
entire region of the flow under consideration the functions v, p, p, p, a 
are continuous but their partial derivatives with respect to x and y 
(all or only some) exhibit jump discontinuities, at least on the curve of 
weak discontinuity OC. The same curve OC appears as a characteristic of 
the second family since the hydrodynamical elements of the stream are con- 
stant. On this line the following relationships will hold in the entire 
region containing the stream 



( 21 ) 
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0 = 8 0 


( 22 ) 


where Qq denotes a quantity defined as 


00 = 



(23) 


From equations (13) , (21), and (22) we easily obtain 


a - 0 


(2k) 


i.e., the flow under consideration is irrotational. By virtue of rela- 
tion (24) the right-hand side of equations ( 15 ) and ( 17 ) vanish and these 
equations can be integrated- As a result of integration of equation (15) 
we obtain the relationship 

3 + cp(v) = constant ( 25 ) 

satisfied along any characteristic of the first family, and as a result 
of integrating equation ( 17 ) we have 


(3 - cp(v) = constant 


( 26 ) 


satisfied along a characteristic of the second family. cp(v) denotes a 
function defined as 


cp(v) 


k + 1 . /k - 1 

arc tani/ 

- 1 / k + 1 



arc tan 



(27) 


Since on the curve of weak discontinuity 0C the quantities v and 3 have 
the values w and 0, respectively, the following relation is satisfied 
along every characteristic of the first family intersecting this line 
and consequently in the entire region of the disturbed stream: 
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P + cp(v) = cp(w) 


( 28 ) 


From equations l6, 2 6, and 28 it immediately follows that the char- 
acteristics of the second family (the curve OC being among these) are 
straight lines since along each of these characteristics the hydrodynam- 
ical elements are constant. 

Making use of these circumstances it is not difficult with the aid 
of equations 28, 26, 22, 21, l6, 10, J, and 5 to construct expressions 
for the functions v, (3, p, p, a in the region of the disturbed flow. 
However, the construction of these expressions is not of great interest 
since our chief interest is centered on the construction of an expres- 
sion for the pressure on the contour K which may be accomplished without 
the use of these expressions for the hydrodynamical elements of the flow. 
Actually equation 28 allows us to determine the velocity v as a function 
of the angle of inclination of this velocity with the x-axis at every 
point of the region filled by the disturbed flow. Since by virtue of 
equation 5 "the angle of inclination of the velocity with respect to the 
x-axis is a given function of x on the contour K there is the possibility 
of using equations 22, 21, 10, 9, and 7 to determine the pressure p as a 
function of x on the flow around a contour. If we limit ourselves to 
the consideration of slightly disturbed flows, i.e., flows whose hydro - 
dynamical elements differ but little from the hydrodynamical elements of 
the undisturbed flow, the expression for the pressure on the flow around 
a contour K may be written in the form of a series. This series has the 
form 


p = p 0 + q a2P k (x) + a 2 p k 2 (x) + a^p^x) + a^p k 4 (x) +.. . . 




( 29 ) 
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where 


q. = 


P 0 w 2 P 0 ^ £ 


a x = 2(M 2 - l)" 1 / 2 


a 2 = (m 2 _ l) -2 ^2 - 2M 2 + ~-i 


a 3 = (M 2 - 1) 


•7/2 


i _ 2A 2 + S(k + 1)M^ + ~ S. + . . 2k2 m 6 + — M 8 

3 3 ^ 8 


at - ( M 2 - !)■ -5(1 - £ m 2 ♦ L+U& (<4 + - 21 - -t 1 ^ 2 M 6 + 

v 3 3 6 12 

lg + 20k - 8k2 + 3k3 m q + -21 - 20k + 3k 2 + 2k? m10 + 3 + 2k - k 2 m i 2 j 

12 48 48 / 


Case IX .- The flow around, a contour K Is situated in such a manner 
that the following inequality is satisfied 


p k (0) > 0 


(30) 


In this case, as is well known, a line of shock discontinuity 0D appears 
(fig- 3) proceeding from the origin 0 and dividing the entire flow under 
consideration in two parts. On one side of this line Is the region of 
the undisturbed stream, and on the other the region in which the fluid 
flows around the fixed frictionless wall. Just as in case I we call the 
flow in the region in which the stream around the fixed frictionless wall 
is accomplished the disturbed flow. In the present case, in contrast to 
case I, the functions v, p, p, p, a exhibit jump discontinuities on the 
shock- line OD. 
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In the region of the disturbed flow these functions must satisfy 
not only equations 2, 5 > and 7 but also the dynamical conditions across 
the shock line. Considering the flow to be only s ligh tly disturbed, these 
conditions may be written in the following form 


y 2 . a2 w2 , V 


2 k - 1 2 k - 1 


( 31 ) 


v = w(l + bjJB + bgP 2 + b ^3 + + . . .) 


( 32 ) 


where 


b]_ = -(M 2 - 


1 

if 2 


b 2 = -(M 2 - 1 )- 2 ^| + M ^ 

7 r * 

b 5 = -(M 2 1)“ i,i+ i M 2 + 2(k _ i) M ^ + Ikl-.l 2 ^ + + ,(.k ± if m 8 

'6 2 k 24 32 " 

b 4 . -( M 2 - Xjfi + 1 M 2 + i£L±J9£ ^ + zLzS2L±22i. „6 + 

\24 8 24 24 

5 + 9k - k 2 + k3 m q + -9 - k - 9k 2 + 9k? 2 d \ 

16 48 J 


— — — 1 + + • . . 
e 0 


( 33 ) 
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where 


4 


l 3 = M^M 2 - 1)' 2 

J 12 


, ^.-a ^m 2 

12 



— = e 0 + e-jB + eoP 2 + . . . (3*0 

dx 


where 


1 

e 0 = (M 2 - l)“ 2 


e l = 


- M^(M 2 - 1)" 2 
4- 


Condition 31 shows that, disregarding the presence of jump discon- 
tinuities in the functions v, p, p, p, a, equation 21, just as in case X, 
is valid throughout the entire region filled by the flow under considera- 
tion. However, condition 22 is not, generally speaking, fulfilled in 
the case now under consideration. However, there is the possibility of 
speaking of satisfying this condition approximately. In fact, consider 
equation 33* Its right-hand side does not contain terms in the first 
and second powers of p. Therefore, for slightly disturbed flows, equa- 
tion 22 may be. regarded as approximately satisfied on the line OD and con- 
sequently throughout the entire region filled by the flow under considera- 
tion. From this it follows that in the region of disturbed flow equa- 
tion 24- may be regarded as approximately satisfied, which means that 



12 


NACA TM 13 94 


equations 25 and 2 6 hold on characteristics. For values of p and v 
near 0 and w, respectively, equation 28 may be written in the form 
of a series 

v = v(l + bq’p + b 2 'p 2 + bj'p^ + b^'p^ + . . .) (35) 


where 


1 

b ± ' = -(M 2 - 1)" 2 = b 1 


V = -(M 2 - 1)" 2 ^| + ifrj = b 2 


7 

b 3 ' = -(M 2 - 1) 2 


bj. ' = -(M 2 - l)~ 5 / r — + — M 2 + T-l .. 2 2 k . M 4 + 2 . — -- 9k yfi + 

4 124 8 24 24 


3 - 2k - 5k 2 + 4 k 3 w8 . -3 + 8k - 7k 2 + 2k^ wl0 

M + M 

32 96 / 


1 . 1 

— + — 

6 2 


M 2 + 2(ic - i)m^ + — 

4 


jl2L±1m6 


12 


6 


Comparing equations 32 and 35 we see that for slightly disturbed 
streams the first may be substituted for the second with good approxi- 
mation. Consequently, for slightly disturbed flows, equation 28 will be 
approximately satisfied along the line OD. Since, on the other hand, 
along each characteristic of the first family equation 25 is approximately 
satisfied, equation 28 will be approximately satisfied throughout the 
entire region of disturbed flow. The approximate expressions for the 
functions v, p, p, p, a are constituted exactly like the accurate expres- 
sions for these functions in case 1. Substituting the approximate expres- 
sion for the function p in the right-hand side of equation 34, we obtain 



NACA 'EM 1394 


13 


a differential equation of the first degree for the approximate deter- 
mination of the form of the shock line. Summing up our considerations 
we can deduce that the accurate results contained in case I can serve as 
approximate results for case II, and further that expression 29 can serve 
as an approximate expression for the pressure on the flow around a contour 
in case II. These same considerations show that there is no sense in cal- 
culating all terms in this expression. It is sufficient to limit ourselves 
to the first two or three terms. 

From all that has been said about cases I and II one may conclude 
that the form of the contour K may be made up in such a manner that art- 
fully constructed shocks may be caused to appear in the region of flow 
around the fixed frietionless wall. In such cases when we pay attention 
to this phenomenon, the results we have obtained are valid, not for the 
entire region of flow around the fixed frictionless wall, but only for 
that part in the neighborhood of the front side of the flow around a 
contour . The fundamental problem of the present work is the construction 
of approximate expressions for the pressure on the flow around a contour 
in case II, with the calculation of the circulation of the flow occasioned 
by the presence of the shock discontinuity OD. In spite of the fact that 
in the case of the presence of circulation it is impossible to integrate 
equations 15 and 1J, there is the possibility, however, of making up such 
combinations of differentials from equations 14, 15, 16, and 17, adding 
to these equations expressions for differentials of the stream function, 
that with the aid of these combinations it is possible to construct expres- 
sions which we shall integrate. Investigations concerning the preceding 
construction constitute the contents of the following section. 


PART II 


Suppose we have a flow corresponding to case II of the preceeding 
section. Assume that in this flow the hydrodynamical elements in the 
region of the disturbed stream differ infinitely little from the hydro- 
dynamical elements in the region of the undisturbed flow. We revamp 
somewhat our notion of the region of disturbed flow. Shortly before we 
agreed to apply this name to the region bounded by the curvelinear 
triangle made up of the curve OCg (contour K), the shock line QCq, and 

the characteristic of the first family C1C2 emerging from the lowest 

point of the contour K (fig. 4) . Taking into consideration equations 5, 
14, 18, and 34, it is not difficult to conclude that with the assumptions 
made just now relative to the hydrodynamical elements the curvelinear 
triangle OC1C2 differs infinitely little from the isoceles straight-line 

triangle 0 , Cq , C2 t (fig- 5) where the equal sides O'Cq' and Cj_'C 2 ' are 
parallel to characteristics of the second and first families in the 
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undisturbed flov. As for the functions p^(x), p, v, p, p, a we assume 

that they all have the properties of differentiability and continuity to 
as many degrees as may be necessary to insure legitimacy of operations 
which are performed upon them. Moreover, we assume that in the flow 
under consideration the infinitesimal quantities p^(x), p-^'fx), (x), 

p, v - w, p - pQ, p - Pq, a - Sq have the same order of magnitude. 

Taking this last group of infinitesimals as fundamental (having unit 
order of magnitude) we shall agree in what follows to adhere to the fol- 
lowing system of notations appearing in investigations involving infinitely 
small quantities. By e m (m being any positive integer) let us denote 

an infinitesimal whose order of magnitude is not less than m. Clearly 
such a mode of notation does not exclude the possibility of several dif- 
ferent infinitesimals being denoted by the same symbol, and, vice versa. 

The same infinitesimal may be denoted by several different symbols. Thus, 
for example, if an infinitesimal a is denoted by ei^, the infinites- 
imal 2a may also be denoted by eij., and, moreover, the infinitesi- 
mals a and 2a may be denoted by e^, e^. 

On an arbitrary characteristic of the second or first family the 
equation 


di|r = pv(sin pdx - cos pdy) 


(36) 


will be satisfied by virtue of equation 8 throughout the entire region 
filled by the flow. Eliminating dx and dy from equations 14, 15 > 
and 36 and taking into account formulas 13 and 21, we arrive at the 
equation 


d(v cos p) + mgd(v sin p) = $ 2 ^ In 9 (37) 

which is satisfied on any characteristic of the first family. Here $2 
denotes the quantity 

a2 |y2sin p cos p - (v^cos^p - a2)m-^j 

®1 = 

k(k - l)v(v2cos2p - a^) (m 2 cos p - sin p) 


(39) 
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On the other hand, having the integral 25 of the equation 

d(v cos p) + m 2 d(v sin p) = 0 (39) 

it is easy to find an integrating factor Lq of this equation, such that 
after multiplying by Lq it may be written in the form 

d[p + cp(v)] = 0 . (40) 

In order to determine we have the obvious relationship 

Li[d(v cos p) + m^v sin p)J = d|p + cp(v)J (4l) 

from which we obtain without difficulty 

Li( m 2 cos P - sin p)vdp = dp (42) 


consequently 


L 


1 = 


1 

v(mp cos p - sin p) 


(^ 5 ) 


If now we multiply both sides of equation 3T by Lq, this equation takes 
the form 


d 


& 



= Hqd In 0 


where Hq denotes the quantity 


(44) 


(.v^cos^p _ a2)nn - v^sin p cos p 

% = - ± 


k(k - l)v2 


(^ 5 ) 
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We denote by H^g the value of at v = w, (3 = 0. We have 


H-1/-1 — 


10 k(k - 1)M 2 


(M 2 - 1) 


1/2 


(46) 


Equation 44 may be rearranged in the following fashion 


d[f3 + cp(v)] = H 10 d in + (H x - Hi 0 )d In 1- (47) 

e 0 e 0 


Wow choose an arbitrary point S in the region of disturbed flow and lead 
a characteristic of the first family through it. We denote the point of 
intersection of this characteristic with the shock line by A (fig. 6). 
Integrating both sides of equation 47 along the above characteristic from 
point A to point S we obtain 


in + f (Hi - Hio)d in ~ 

8 0 ) J AS e 0 

(48) 

where (3 S , v 8 , 0 S denote, respectively, the values of (3, v, 0 at 
the point S and p a , v a , 0 a denote the values of these quantities at 
the point A. Taking account of equation (32) we have 


Ps + <P(v s ) - P a “ T(v a ) - H 10 (ln Q 




6s 

'0 


v a = w (! + b lPa + b 2Pa 2 + b 3Pa ^ + b 4Pa^ + ■ ■ ■ ) 

(49) 


We introduce the quantity v R -| defined with the help of the expansion 35 
in the following fashion 


v a i = w(l + bxpa + b 2 P a 2 + tj’Pa 3 + b 4'Pa^ + • ■ «) 


(cf eq. (35) -Tr.) 


( 50 ) 
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By this definition of the quantity v a i we have 

Pa + «P(val) = <P(») (51) 

With the help of formulas kS, 50* and 51 we rearrange the expression 
Pa + cp(v a ) in the following manner 

Pa + «p(v a ) = P a + CpCvgo.) + cp(v a ) - 9^^) 

= cp(w) + ?(v a ) - cp(val) 

= cp(v) + cp’(w)[(v a - w) - (v^ - w)] + 

| 9 '"(w) £(v a - w) 2 - (v al - w) 2 J + . . . 

= cp(w) + wcp'(v) jjb 5 - bj'JPa 3 + (ty. - bi^ ' ) P a ^j + 

w^p"(w)b 1 (b3 - bj^p^ + (52) 


Calculating cp * (v ) , 


<p" (v) we obtain 


9 ' (w) = - 


1 

wb^ 


2b o 

9" (w) = — — 


(55) 

(5*0 
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Using formulas 52, 53) and 54 ve easily find 


Pa + <p(v a ) = <P(w) 




2bo, 

_(b 5 - b 3 >) 


— (b4 - b4‘) 
b l 



+ G 


5 


(55) 


Now pass a stream line through the point S and denote by P the point of 
intersection of this line with the shock line. Since the stream function C 
is constant along this line we have 


e B - 


(56) 


where 0p denotes the value of 0 at the point p. Taking logarithms of 
both sides of equation (33) we obtain 


In 



+ Z^'p^ + + 


(57) 


Since the values of the coefficients Zjj.', * , . . . will not be needed 
in what follows, we shall not calculate them. 

Using formulas ( 56 ) and (57) we easily see that 


3* — - In — = Z 3 (p p 5 - Pa 3 ) + V<Pp 4 “ Pa 4 ) + £ 5 (58) 

9 0 e 0 


where p p denotes the value of p at the point p. Assuming that the 

mean value theorem is applicable to the integral arising from the right- 
hand side of equation (48), we easily find^- 


-*-Instead, take a slightly more general assumption admitting the part 
AS of the characteristic under consideration to be divided in the same 
finite number of parts in such manner that on each part the mean value 
theorem can be applied to the integral under investigation. 
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(Hq - Hio)d la i-= 13 ( 5 ! - E 10 )(p p 3 
e 0 


Pa 3 ) + e 5 


(59) 


(in consequence of this equation one must keep in mind that Hq - H 10 = 6 1 
Here Hq denotes the value of Hq at some point on the characteristic 
under consideration between the points A and S. 

Using relations (55 ) > (5&)> and (59) we write equation (48) in the 
following form 


3s + 9(vs) ■ 9( w ) ~ r^ b 5 " b 3')Pa 3 + 

b l 



b V) 


A( bit - b^*) 

b l 



+ 


HlO^Pp 3 - 3a 3 ) + HloV(3p^ - 3a 11 -) + 

^ 3(^1 “ a lo)(3p 3 - 3a 3 ) + e 5 (^°) 


We denote by B the intersection of the characteristic of the first 
family under consideration with the contour K. Applying formula (60) to 
the point B (which is possible, since the point S was chosen arbitrarily) 
we obtain 


3b + <P(v b ) = <p(w) - i-(b 5 - b 5 ')p a 5 + 


2b 


b l‘ 


(b 5 - b 5 ') - -i-(b4 - b] i .') 


3a 14 - + 


%0 2 3(3o 3 - 3a 3 ) + H 1 0 V(3 0 4 - Pa^) + 


l 5 (Hq - Hqo) (p 0 ? - 


(61) 
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where p b , v^, denote, respectively, the values of p and v at the 
point B and po denotes the value of p at the point 0. 

We now proceed to the derivation of an expression for p a . From 
formula (60) we have 


Ps + ‘POO = *P(w) + £5 (62) 

From equation ( 62 ), using formula (35) we obtain 

Vg = w(l + bxps + b 2 p s 2 ) + £3 ~ (63) 

and denoting by m2 B the value of m 2 at the point S we obtain, by 
using formulas (19) and (63) 

m 2s = (M 2 ~ 1) _1/ ^ 2 + — -g - M^(M 2 - l)" 2 p s + e 2 


= e 0 + 2eqP s +' e 2 


(64) 


Analagous to the derivation of equation (47) , which holds on character- 
istics of the first family) we may derive equation 


djp - cp(v)] = H^d In ~ + (H2 - In £- 


8 0 


(65) 


which is valid on characteristics of the second family. Here S 2 denotes 
the function defined as 


H 2 = 


(v2cos2p - a 2 )!^ - v 2 sin p cos p 


k(k - l) v 2 


( 66 ) 


and H2 q denotes the value of this function at p = 0, v = w. 
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Now pass a characteristic of the second family through the point S 
and denote by Q its intersection with the contour K. Integrating both 
sides of equation (65) along this characteristic from the point Q, to the 
point S we obtain 


Ps “ Pa " 


19 (vs) 


- 9(v q ) 




0 , 


= Hao In — - In + 


V 


So 




QS 


(H2 - H20)<i In ~ 

°0 


( 67 ) 


where p q , v q , 0 q denote respective the values of p, v, 0 at the 
point Q. Since the contour K is a stream line we have 

0 q = e(°) (68) 


where 0 ^°) denotes the value bf 0 at the point 0 . Assuming that the 
mean value theorem can be applied to the integral arising from the right- 
hand side of equation (67) we easily find, with the aid of formulas (56), 
(57), and (68) 


Ps “ Pa “ 


9(v s ) - 9(vq) 


“ e 3 


(69) 


Applying formula (62) at the point Q, we have 


Pq + 9(v q ) = 9 00 + 


( 70 ) 


Eliminating cp(w) from equations ( 62 ) and ( 70 ) we arrive at the fol- 
lowing equation 


Ps - Pq + 



(71) 


9(v s ) - q>(v q ) 
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Families (69) and (71) give 


P s = Pq + ( 72 ) 

On the shock line we take an arbitrary point F (fig. 7) and paBS 
through it a characteristic of the second family in the region of the 
disturbed flow and we denote by p£, respectively, the values p 

and m2 at the point F. Applying formula (64) at the point F we obtain 

m 2f = e 0 + 2e iPf + e 2 ( 75 ) 

We denote by ( — ) the slope of the tangent to the shock line at the 
Vdx/f 

point F. From equations (34) we have 


dx f 


= e 0 + e lPf + e 2 


(74) 


Comparing formulas (73) and (74) we see that the characteristic of the 
second family passing through F and the shock line at this intersection 
make an infinitesimal angle with each other moreover, if 


Pf > 0 


(75) 


the slope of the characteristic of the second family is greater than the 
slope of the shock line at the point F. 2 


^It is easy to show that if the shock line is unbroken and moreover 
condition (30) is satisfied the inequality p < 0 is impossible on this 
line. As a matter of fact, in the opposite case the shock line is broken 
since with p < 0 condition (34) must be replaced by the following con- 
dition in virtue of Tsemplen's theorem 


ay 

dx 


- ( e 0 " e lP + e 2P 2 ~ 


■ • •) 
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Denoting by L the intersection of the characteristic of the second 
family under consideration with the contour K and by x x abscissa of this 

point, we have 


*1 = e i (76) 

Let p 2 , denote the value of p at the point L. Using equations ( 5 ) 
and ( 76 ) and MacLauren's formula we obtain 

P x = 0 k (O) + p k '(0)x x + e 3 (77) 

Applying formula ( 72 ) at the point F we obtain 

Pf = P x + (78) 

As a consequence of equations ( 77 ) and ( 78 ) 

= P k (0) + P k l (0)x z + (79) 


Since the point F was chosen arbitrarily on the shock line by use of 
equations (74) , ( 76 ) , and ( 79 ) we can obtain the following differential 
equation for the shock line 


jjjf = e o + e iM°) 


+ e 2 


( 80 ) 


Consequently the equation of the shock line may be written in the fol- 
lowing form 


y = 



+ e lPk(°) 


X + €£ 


( 81 ) 
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Applying formulas (64) and (72) to an arbitrary point situated on the 
characteristic of the second family PL we easily obtain the differential 
equation of this line from the following form 


+ 2e-j_3£ t £ 2 ( 82 ) 


this equation may be written 


— = e 0 + 2623^(0} + £2 (83) 

dx 


Consequently the equation of the characteristic PL may be written in the 
form 

y = y 1 + Je 0 + (x - + e 2 (84) 


where y^ denotes the ordinate of the point L. 


On the other hand, taking account of formulas (3) and ( 76 ) we have 


<47 

— = e° 
dx 


Employing formulas ( 76 ) and (77) 


Tl = / tan 3 k (x)dx = e 2 

J 0 


(85) 


Employing formulas ( 85 ) and ( 76 ) we may write equation (84) in the form 


y = 


e 0 + 2e iPk(°) 


x “ e o x l + e 2 


( 86 ) 


Applying formulas (8l) and (86) at the point F and denoting by Xf, y f 
the coordinates of this point we obtain 
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y f = e 0 + e 1 p l5 .(0) 


Xf + e 2 


(87) 


y f = e o + 2 e iPk(°) 


Xf - e 0 x x + e 2 


From equation (87) we easily obtain 


xi = — x f p k (0) + e 2 


( 88 ) 


Replacing x^ in the right-hand side of equation (79) by the expression 
in formula (88) we obtain 


Pf = Pk(®^ + x f‘P k (0)P] C I ( 0 ) + e, 
e 0 


(89) 


We denote by x a , y a the coordinates of the point A and by x^, y^ 

the coordinates of the point B. Applying formula (89) at the point A 
we arrive at the following result 


p a = p k (0) + — XaPkCOPfc'CO) + s 3 
e 0 


(90) 


We now express x a in terms of x^. To this end, using formulas (l4) 

and (l8) , we write the differential equation for characteristics of the 
second family in the following fashion 
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Employing formula (91) we write the equation for the characteristic AB 
of the first family in the form 


y = Yb - e o( x - xb) + e ± 


(92) 


Taking account of formula (3) we have 


yt - 



tan 0^(x)dx = 


e l 


(93) 


consequently equation (92) may be written 

y = -e 0 (x - x b ) + e x (9*0 

On the other hand, equation (8l) for the 6hock line may be written in 
the form 


y = eo x + 

and applying formulas (9*0 and (95) at the point A we obtain 

\ 

y a = " e o( x a - x b) + .61 

> 

y a = e o x a + £ 1 

J 

From equation (96) we easily find 


(95) 


(96) 


x 


a 



e l 


(97) 
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And consequently 


P a = Pk(0) + “=- ^bPkCOiPk'CO) + e 3 

2e 0 


(98) 


Substituting this expression for p a in the right-hand side of equa- 
tion (6l) and substituting (Sq for p k (o) in the fundamental formula ( 5 ) 
we obtain 


Ob + <P(vb) = «P(w) - — (b 5 - b 5 ')j3 k 5(o) - 


2b c 


— (bl). - blj. 1 ) “( b 3 ~ b 3’) 


Pk (°) 


3e l 

(b 5 - b 3 ' 

2e 0 

1 

cr 

H 


+ h 10 i 5 


x bPk 3 (°)Pk'(°) + e 5 (99) 


Employing relation (35), we easily obtain from equation (99) 

/ 

v b = w • 1 + bqp b + b£0b 2 + b 3'0b^ + h^'Pb^ + ( b 3 " b 3*)Pk^( 0 ) + 


2b r 


b 4 - bit* (bj - bj 1 ) 


0 k 4 (o) + — ( b 3 - b 3 , )0 k 3 (o)0 b + 


— i[b 5 ■- bj* + H lo ^3b ; Ljx b 0k 5 ( o )Pk , (°) 


( + € c 


(100) 
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Substituting v^, p^, *b for v > p k (x) and. x, respectively, in 

formula (100) we arrive at the following final expression for the velocity 
on the contour K: 


v = w Jl + b3JB k (x) + b 2 P k 2 (x) + b3'p k 3( x ) + ty/p^Cx) + (bj - b 5 ')p k 5 (0) + 


2b, 


t>4 - b4* - — ( b 3 - b 3') 


P k 4 (°) + “( b 3 - V)Pk 3 (°)3 k (x) + 


3e 


2e 0 


~[ b 5 - b 3* + Hi 0 i 3 bjxp k 3 (0)p k '(0) 


+ €c 


(101) 


We now proceed to the derivation of formulas from which the pressure 
on the contour K can be calculated. Clearly 



( 102 ) 


and moreover, on the contour K the following equation holds 


JL m 0 (o) 

p k 


(103) 


Employing formulas (7), (10), (21), (23), (102), and (103) we easily obtain 
the following expression for the pressure on the contour K 


1 


e(°) 

k-l - 

, 00 _ 



k - 1 



P = Po. 


2 


k 

I k-1 


(104) 
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On the other hand, by virtue of equations ( 5 ) and (33) the following 
equation holds 


= 1 + 1 ^( 0 ) + 11 ^( 0 ) + . . . ( 105 ) 

0 O 


Substituting the expressions for v and obtained in formulas (lOl) 

e 0 

and (105 ) ) respectively, in the right-hand side of equation (104) we obtain 
after elementary transformations the desired formula for the calculation 
of the pressure on the contour K 


P - Po + q.[a-iP k (x) + a23 k 2 (x) + a^p^x) + a^^x) + a ld p k 5 (0) + 

a 2d0k^(°) + a 3dPk 5 ^°^k^ + a 4a.P k 5 (°)Pk' C 0 ) 3 ^ + e 5 ( lo6 ) 


where 


a ld = -2(b 5 - bj') - 


21 , 


k(k - i)m 2 


k + 1 


M^M 2 - 1)' i(- | + M 2 + ^ 


a24 = -2(blj. - bV) + — ^(b 3 - b 5 ') - 


21, 


k(k - 1)M C 


= M^(M 2 - i)-5/l k+J: + 2 - 1 - 2L I m 2 + - 10 - + - 6k2 r m 1 ^ + 

\ 2 k 8 


9 - 7 k 2 + 2 k? ^ + -3 + k + 3 k 2 - k? m q’ 


16 


32 
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a 3d = ( b 3 - b 3*) j^ 2b l - 2b x - —^ ) + 


4bg\ 2l^bj_ 


bl/ k - 1 


= 2^(M 2 - l)“ 5 /- 1 + L±.. 2k . ~_1 I E M 2 + 

6 24 


-4 + 3k + 6k 2 - k3 ^ _ 3 - 7k - 7k 2 + 3k3 


24 


96 


M 6 


3e 


a 4d = - ( b 3 “ b 3‘ + H l0*3 b l) 


= — +- 1 ) M 8 (M 2 - l)~^(-l + ? ■ ~ k M 2 + - k 5 mA 

16 y 2 8 ) 


Por x = 0 the formula (106) takes the form 

P = P 0 + a[»iPk(0) + a 2Pk 2 ^ 0 ^ + a 3'Pk 5 ^ 0 ^ + a 4*Pk^ 0 )] + e 3 (iOT) 


where a^ 1 = a^ + a^, a^' = a^ + ag^ + a^. 

Formula (107) may be used for the calculation of the pressure on a 
flat plate which Is inclined at an angle p^(0) to the undisturbed flow. 

In order to single out of the right-hand side of equation (106) those 
terms which depend exclusively on the presence of the shock in front of 
the contour K, we add to the contour K under consideration an arc 0*0 
of finite length in such a manner that this arc is tangent to K at the 
point 0 and is parallel to the x-axis at O' (fig. 8). Since the flow 
around such an additional contour is accomplished without the appearance 
of shocks (we suppose that the angle between the direction of flow and the 
x-axis and the derivative of this angle with respect to x are both 
infinitely small ) , formula (29) may be employed in the calculation of the 
pressure on this contour. Comparing formulas (29) and (106) and denoting 
by dp s toss tbe pressure resulting from the presence of the shock front, 
we obtain 
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^?stoss 


9 j a ldPk'^°) + a 2dPk^(°) + a 3dP^(°)Pi 5; ( x ) + 
a 4dPk 5 ^°^k , ^°^ x ] + e 5 


(108) 


We may, in. turn single out of the expression for ^?> B toss "fc erm 

depending solely on the vorticity caused by the presence of the shock. 

In order to do this -we add to the contour K(0C 2 ) a straight-line segment 

tangent at the point 0 (segment O'O in fig. 9)* With the contour 0*OC 2 
a shock is formed, but the shock line O'C^ is straight so that vortex 
formation is absent. Calculating the pressure on the portion 0C2 of the 
contour 0 * 002 * ' we obtain 


P = P 0 + q[a 1 p k (x) + a 2 p k 2 (x) + a^p^x) + a^p^x) + 

a ldPk^°) + a 2dPk^^ + a 3d^k^(^)Pk^ x ^| + e 5 (l°9) 


Comparing formulas (106) and (109) and denoting by ^P ro ^ "the pressure 
due to vortex formation caused by the shock, we obtain 

Z^rot = ( l a 4dPk 5 ( 0 )Pk , ( 0 ) x + e 5 (H°) 


PARI III 


We now apply the results obtained to the calculation of the lifting 
force and head resistance of a flat wing with sharp front and rear edges 
placed in a supersonic stream having constant hydrodynamical elements. 

We place the origin 0 at the front edge of the wing and arrange the 
coordinate system so that the positive x-axis corresponds to the direction 
of the velocity of the undisturbed flow and measure angles in the manner 
used heretofore. Segment OC 2 connecting the front and rear edges (fig. 10) 
will be called the chord of the wing as in the theory of wihgs. The length 
of this curve will be denoted by I and the angle it makes with the x-axis 
by p. 
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The form of the wings we are investigating is defined by a pair of 
contours like that investigated in the preceding section, possessing a 
pair of common points 0, C2* Comparing ordinates of points on these 
contours having the same abscissa, we call the upper contour K u that 

contour of which every point on the ordinate is greater than the corre- 
sponding point on the ordinate of the other contour, which we call the 
lower contour Kj . The function pj^(x) for the upper contour we denote 
by Picu(x) and for the lower by Pkj(x). 

We choose an arbitrary point A on the chord of the wing and denote 
the distance OA by t. Through A we pass a straight line perpendicular 
to the chord of the wing and denote by A^ and A j, respectively, the 

intersections of this straight line with the upper and lower contours. 

With the point A u we associate a unit tangent vector tu and at the point Aj 

a unit tangent vector t j . The vectors t^ and tj will be directed in such 

a manner that their projections on the direction 0 C 2 are positive. We 
denote by p u and p^, respectively, the angles these vectors make with 
the vector OC 2 . Clearly p u , p^ may be regarded as functions of t. 

We denote by p u0 and p J0 the values of p u and at the point 0 

and the values of the derivatives of p u and Pj with respect to t 
at the point 0 by p u0 ' and P^q** 

The abscissa x of A u may be calculated from the formula 


t n t 

x = t cos p - Bin p / tan P u dt 

J 0 


( 111 ) 


and that of the point Aj from the formula 


x = t cos p - sin p / tan P^dt 


L 


(112) 


The value of the functions Pj^Cx) at the point A,^ is determined by the 
relation 


(H3) 


Pku( x ) = P + Pu 
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and. the value of the function )3 tj.^(x) 


at the point A-i by the relation 


Pk l( x ) = P + Pi 


(114) 


Moreover we have the relations 



tan (3 u cLt = 0 


(115) 



tan fSjdt = 0 


( 116 ) 


We assume that j3 and also (3 U and and other derivatives with 

respect to t are infinitesimal quantities. From equations (115) and 
( 116 ) we easily obtain 



(117) 


( 118 ) 


Proceeding now to the calculation of the lifting force and head resistance 
of the wing under consideration, we remark that on the top side of the 
wing a shock appears when and only when 


P + P u 0 > 0 


(119) 


and at the bottom side when and only when 


P + Pio < 0 


( 120 ) 
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We introduce the quantities aq u , 8 -2m> a 3u> a 4u ^- e ^i ne( i as follows 


a lu = a ld I 


a 2u = a 2d 
a 3u = a 3d 


if 


p + > 0 


> 


a 4u = & kd J 


a lu - a 2u = a 3u “ a ku = 0 if P + Ho * 0 


( 121 ) 


( 122 ) 


In an analogous way we define' the quantities a-^, & 2 .i> a 3 1 > a ^j 


a 12 
a 2 1 
a 3 1 
a kl 


= a 


a ld 

2d 

a 3d 

a 4d 


if 0 + Pi 0 < 0 


(123) 


a ll ” a 22 = a 32 ~ a lj.j - 0 P + PjO - ® (124) 


Denoting by p u the pressure on the upper contour and by pj 
the pressure on the lower contour Kq we easily obtain, with the help of 
formulas (106), (121), (122), (123), (124) 

P u = P 0 + QL (aiPfaik) + a 2Pku 2 ^ x ) + a 5Pku 3 ^ x ) + a 4i 3 ku it ( x ) + 

a luPku 3 (°) + a 2uPkuV°) + a 3uPku 3 ^°)Pku( x ) + a 4uPku 3 ^°)Pku' (0)x ] + e 5 

(125) 
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pz = po + q-^iPkzCx) + a 2Pki 2 ^ “ a 3Pkz 5 ( x ) + Vk z^ x ) - 

a lZPkl 5 (Q) + a 2ZPkl^°) + a 3^ P kZ 3 ^°^kl^ + a 4z P kZ 5 ( 0 ^kl'( 0 ) x ] + e 5 


(126) 

— ^ 

Let P denote the resultant vector of the hydrodynamical force 
acting on a unit length of the wing under consideration. We have then 


P = 



(127) 


where n denotes a unit vector normal to the contour of the wing and 
directed inwards. 


We introduce the dimensionless coefficient of the lifting force Cy- 
an! the dimensionless coefficient of head resistance C x . These coeffi- 
cients are defined by the formulas 



(128) 


C x = 



(129) 


where Py and P x denote the projections of the vector P on the 
x- and y-axes, respectively. Prom formulas ( 127 ), (128), (129) "we have 


Cy = - — [ Pudx + — f Pjdx (130) 

y qT Jo qT J 0 
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c x = + cp J'o Pu tSn Pku(x)dX “ ^ X Pl tan Pkl(x)dX (151) 

where 1 denotes the abscissa of the point C 2 . With the aid of formu- 
las (111), (112), (113), (2JJ0, (H-7), ( 118 ) , (125), ( 126 ), (W°)i 
nri ri (131) after a few elementary transformations we obtain 


Gy = Cyl + Cy 2 + Cy 3 + C y 4 + €5 


(132) 


where 


Gyi = -2a 1 (3 

V - - T / 0 T (*»* - ^1“ 


3-3 
C y5 = (ai - 2a3)p^ - a i u (P + Pud) " a lZ^ + Pzo) + 


"y3 


i(al - 3a,)? (Pu" ♦ P» 8 )* ♦ - - 3 ) / 0 (P“ 5 + 

c y4 * -a^S + + a 2l($ + Plo)** " a 3 u P(P + Put)) 3 + 

ajtfCS + PlO> 3 - 1 «m<P + 3 u0 ) 3 PaD , '+ 1 a 4j(P + PlO ) 3 lW + 

i^| aa - 6a^p 2 u 2 - Pl 2 )<*t + ^ a 2 “ ka kj P ~ 


Pi 5 ) 4t 
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Ox - &x2 + c x3 + O x ij. + e 5 


(133) 


0x2 = 2a x p 2 + y / ( p u 2 + Pz 2 ) dt 


Cx3 * f Q (Pu 2 - Pi j" ♦ £ / o (p u 3 - Pi 3 )dt 

0x4 = ^2a^ - + a lu^P + Puo)^ + a iiP(P + PIO) 3 + 

l( 6a 3 - t)p 2 f* (^ 2 + pi 2 )" + i(*“3 - y)? / 0 T (pu 3 + Pi 

?( a 3 + t) XV + p V 


Let us consider a numerical example. 


k = 1.405, M = 1-5 , Pq = 1*033 kg/cm 2 

q - P °^ = I.633 kg/cm 2 

a x = 2(M 2 - 1 r 1 / 2 = I.789 
a 2 = (M 2 - 1)“ 2 (2 - 3 A 2 - + 1.203M**-) = 2.296 
a 5 = (m 2 - 1) “7/ 2 (i. 333 - a^ 2 + 4.00 - 1.81346 + 
o.4ooai 8 ) = 3*082 

a^ = (M 2 - 1) (O.6667 - O.6667M 2 + 5.61^ - 3.82^ + 
2.96SM 8 - 0.7840M K) + 0 . 0799 2M 12 ) = 8.290 


(Equations continued on next page). 



38 


NACA TM 1 39k 


aid. = 1*203M^(M 2 - l)" 7//2 ( -0.3333 + 0 . 265^2 _ o. 03271 M^) = 0 . 2? 66 ] 

a 2d = M^(M 2 - 1)" 5 (-1.203 + 1-317M 2 - + O.Ol^M 6 + 

o.o 485^ 8 ) = o.Wt-8 

a 3d ” ^(M 2 - 1) “^(-0.1(008 - 0 . 00234 2 + 0.3869^ - 0.1285^) 
ai|d. = O^ISM^M 2 - 1)“ 5 (-1 + 0.79754 2 - 0.09813^) = 0.9035 

Let us take as the functions p^ p^ 


135 \ 

^ I cone .1 


= 0.3318 


P u = -2(3 + 



H - 0 


( 136 ) 


Moreover we assume that 


P < 0 


(137) 


The form and position of the profile of the wing, determined by equa- 
tions ( 138 ) and condition (l37)> is shown in figure 11. It is easily 
seen that the straight line S-j_S 2 drawn perpendicular to the wing through 

its mid point is the axis of symmetry of the profile under consideration. 
From equations ( 136 ) and condition (137) we have 


3 + PuO = -P > 0 

► 

P + p I0 = P< 0 


(138) 
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Consequently, 


a lu = a lZ = a id ■ 0.2766 


a 2u = a 2l - a 2d = O.Wl-8 


l 5u = a 3Z = a 3d “ 0, 33l8 


a4 u = a^i = al^_ = 0.9035 


Using formulas (132), (133 ) > (l35)> (136), and (139) we obtain 


Cy = -2a;jj - - a 2 p 2 + a x - 6a^\ 0 3 + 


10 36 

_ a 2 _ __ alf + 2a 5d + + £5 


-3.578(3 - 3.06l| 2 - 14.32P 5 - 82.73^ + 


°x = ~ a lP 2 + + + ^ a 3J ^ + e 5 


5.963P 2 + 9.184-p 3 + 4o.8p^ + . . . 


0 =--T 

36 


T = 100 cm 
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Then 


G y = 0.2936 
Cx = 0.04l68 


(142) 


Knowing Cy, Cx, T, and q, we easily obtain 

Py = qTCy = 47*9 kg/cm 




P x = qTCx = 6.8l kg/cm 


(143) 


SUMMARY 


In the present work the problem of a flow of stream of ideal gas 
around a thin wing at small angles of attack is investigated, this stream 
being supposed to be two-dimensional, stationary, supersonic and deprived 
of heat-conduction. 

In the initial part of the work, the problem is stated, and the well- 
known results obtained by Ackeret, Prandtl, and Busemann are cited. These 
results, as known, are obtained on the basis of the potential supersonic 
streams theory, which is founded on the existence of integrable combina- 
tions of characteristics of differential equations concerning this problem, 
and in which some peculiarities of the dynamical conditions on the line of 
the shock wave are utilized. 

In the second part the approximate solution of the problem is given 
with an allowance for vortex- format ion caused by the change of entropy 
along the shock wave, when receding from the leading edge of the wing, 
near which this shock wave is formed. For this purpose differential 
equations of characteristics non admitting integrable combinations are 
to be dealt with. The solution is obtained by means of a special method, 
which enables us to find the approximate integrable combinations of dif- 
ferential equations of the characteristics. The obtained combinations 
let us receive the approximate formula of pressure in any point of the 
contour of the wing investigated. From this formula the term is easily 
segregated depending exclusively on the vortex formulation, caused by 
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the change of entropy along the shock wave. The characteristic dis- 
tinction of this term of the obtained formula of pressure from the other 
ones, is that it includes the curvature of the wing contour at the leading 
edge and the distance from this edge up to the element of the wing for 
which the pressure is calculated. 

In the third part of the work the expressions for lift and drag 
coefficients of the wing are given, on the base of the formula of pres- 
sure obtained above. In conclusion a numerical example is studied. 


Translated by R. Shaw 

Institute of Mathematical Sciences 
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Figure 5 . 



Figure 6 . 








